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The process of a light neutral vector meson electroproduction is studied in the framework of QCD factorization 
in which the amplitude factorizes in a convolution of the nonperturbative meson distribution amplitude and the 
generalized parton densities with the perturbatively calculable hard-scattering amplitudes. We derive a complete 
set of hard-scattering amplitudes at next-to-leading order (NLO) for the production of vector mesons, V = p°, lj, <j>. 



1. The process of elastic neutral vector meson 
electroproduction on a nucleoli, 



7*(q)N(p)^V(q')N(p'), 



(1) 



where V — p°, ui, 4>, was studied in fix target and 
in HERA collider experiments. The primary mo- 
tivation for the strong interest in this process (and 
in the similar process of heavy quarkonium pro- 
duction) is that it can potentially serve to con- 
strain the gluon density in a nucleon [1,2]. On 
the theoretical side, the large negative virtuality 
of the photon, q 2 = —Q 2 , provides a hard scale 
for the process which justifies the application of 
QCD factorization methods that allow to sepa- 
rate the contributions to the amplitude coming 
from different scales. The factorization theorem 
[3] states that in a scaling limit, Q 2 — ► oo and 
%Bj = Q 2 /2(p • q) fixed, a vector meson is pro- 
duced in the longitudinally polarized state by the 
longitudinally polarized photon and that the am- 
plitude of the process (1) is given by a convo- 
lution of the nonperturbative meson distribution 
amplitude (DA) and the generalized parton den- 
sities (GPDs) with the perturbatively calculable 
hard-scattering amplitudes. In this contribution 
we present the results of our calculation of the 
hard-scattering amplitudes at NLO. 
2. p 2 = p' 2 = m 2 N and q' 2 = m 2 Ml where rtiN 



and tum are a proton mass and a meson mass 
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respectively. The invariant cm. energy s 7 * p = 
[q+pf = W 2 . We define 
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P + p' 
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l M ' x Bj 



t = A 2 . 
Q 2 



W 2 + Q 2 

We introduce two light-cone vectors 
n 2 + = n 2 _ = , n + n_ = 1 . 
Any vector a is decomposed as 



(2) 



(3) 



a M = a+n^ + a~nt + a± , a 2 = 2 a+oT - a 2 . (4) 

We choose the light-cone vectors in such a way 
that 



p=(l + £)Wn+ + 



P 



' - (1 - 0Wn + + 



2(1 + £)W 
A 2 ) 



(5) 



2(1 -g)W 



We are interested in the kinematic region where 
the invariant transferred momentum, t, is small, 
much smaller than Q 2 . In the scaling limit vari- 
able £ which parametrizes the plus component of 
the momentum transfer equals £ = XBj / (2—XBj)- 
GPDs are defined as the matrix element of the 
light-cone quark and gluon operators [4] : 
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H q (x^,t)u(p') ji_u(p) 

(6) 
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In both cases the insertion of the path-ordered 
gauge factor between the field operators is im- 
plied. Momentum fraction x, — 1 < x < 1, 
parametrizes parton momenta with respect to the 
symmetric momentum P = (p+p')/2. In the for- 
ward limit, p' = p, the contributions proportional 
to the functions £ q (x,^,t) and £ 9 (x,£,i) vanish, 
the distributions Ti q (x, £, t) and H 9 (x, £, t) reduce 
to the ordinary quark and gluon densities: 



H q (x,0,0) = q(x) for x > 0, 
H q (x,0,0) = -q(-x) for x<0; 
H 9 (x,0,0) = xg(x) for x > 0. 



(8) 



Note that gluon GPD is even function of x, 
W{x,^t) = W{-x,i,t). 

The meson DA <fiv(z) is defined by the follow- 
ing relation 

(o\q(yhM-y)\VL(p)) y ^ (9) 
i 

= p^f v Jdze^- 1 ^4> v (z). 



It is normalized to unity / <p v (z)dz = 1 . Here z 

o 

is a light-cone fraction of a quark, /y is a meson 
dimensional coupling constant known from V — > 



decay, in particular, /„ 



27rV47ra /v 



198 ± 7MeV. 
i i 

dz <pv{z) J dx 



Q v (r g (z, x) F"(x, & t) + T (+) (z, x)FW ( X; t)) 



+ X)ej:T,(« )I )F(+)( I) f 1 t) 



(10) 



Here the dependence of the GPDs, DA and the 
hard-scattering amplitudes on factorization scale 
jj,F is suppressed for shortness. Since we consider 
the leading hclicity non-flip amplitude, in eq. (10) 
the hard-scattering amplitudes do not depend on 
t. The account of this dependence would lead to 
the power suppressed, ~ t/Q, contribution. 2 a is 
a fine structure constant, N c = 3 is the number 
of QCD colors. Qv depends on the meson fla- 
vor content. If one assumes it is -j=(\uu) — \dd)), 
-j^(\uu) + \dd)) and \ss) for p, ui and <f> respec- 
tively, than Q p = -±= , Q w = -±= and Qj, = =±, 
the sum in the last term of (10) runs over q = u,d 
for V = p, ui and q — s for V — <ft and 
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F q( -+\x,£,t) = F q (x,£,t) - F q (-x,£,t) de- 
notes a singlet quark GPD, F( + \x,£,t) = 
T, Q =uds F9{+) ( x '^ t ) stands for the sum of all 
light flavors. 

Due to odd C- parity of a vector meson tfiv (z) = 
</>y(l — z). Moreover, since V and 7* have 
the same C- parities, 7* — > V transition se- 
lects the C-even gluon and singlet quark contri- 
butions, whereas the C-odd quark combination 
F q (-\x,£,t) = F q (x,£,t) + F q (-x,£,t) decou- 
ples in (10). 

3. Below we present the results of our calcula- 
tion of the hard-scattering amplitudes in the MS 
scheme. 

T q (z, x) may be obtained by the following sub- 
stitution from the known results for a pion EM 
formfactor, see also [81, 



T q (z,x) = {T(z,^A- l , 



2£ 



2 For the analysis of other helicity amplitudes see [5-7]. 
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+6 ln(t)ti) 



ln(w) ln(u) 28 



/5 /Q 2 
+/3 - ln(^) - In 

+c2 / + y _ L , 

\ (U — m) ,:1 

+Lz 2 (t>) — Li 2 {u) + ln(w) ln(w) — ln(u) ln(u)] 

„ (v + u — 2vu) In vu n . , , . , 

+2^ ? + 2 Lz 2 (u) + Li 2 (u) 

[v — uy 

-Li 2 (u) - Li 2 (v) + ln(v) Ln(it) + ln(u) ln(v)} 

vu\n(vu) a , ,_. , ,_. 20 
+4— L 5 i_41n(«)ln(«)-y 

Here and below we use a shorthand notation 
u = 1 — u for any light-cone fraction, fin is a 
renormalization scale for a strong coupling, f% = 
iig^ — ^21, Uf is the effective number of quark 

flavors. Cf = Li 2 («) = - J f ln (l ~ l )- 

o 

Also we denote 
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T( + )(z,a;) starts from NLO 



(8tt)zz 



here 
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-(y + l)ln(y + l))]-^ 

/J y + z 
yln(-y) + (y + l)ln(y + l) 

y{y + i) 

y(y+l) + (y + z) 2 
(y + z) 2 

{*-»*}' 



(15) 



where we introduced two auxiliary functions 

R{z,y) = zln(-y)+zln{y+l) 

+zln(z) + zln(z) , (16) 

H(z, y) = Li 2 {y + 1) - Li 2 (-y) + Li 2 {z) 
-Li 2 {z) + ln(-y) ln(z) - ln(y + 1) ln(z) . (17) 

For the gluonic contribution we obtain 



Tg(Z,X) 



: (a; + £ - ie)(x -£ + ie) 



+ 4ir 9 {' 2£ 



le 



(18) 



Q' 2 



where 

X s( z,y) = {(ln^-]-l 

A) 2( Cl -c 2 )(y 2 + (y + l) 2 ) / 
2 y(y + l) I 

(y + l)ln(y + l)-yln(-y)) 

+ ci / ym(-y) + (y + 1) ln(y + 1) 



2 V y + i 



+ci (^- + 2zln(z) 

2 V Uh/ / 2(y + z) 

(3 Cl -4c 2 ) /ym 2 (-y) 
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4 V 2/ + 1 
(y + l)ln 2 (y + l) 



(ln(-y) + ln(y + 1) 
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Ci ^zln(z) - + 2c 2 

-(ci -c 2 )(ln(zz) -2) 
(y + l)ln(y + l) 



+ ci(l + 3z)ln(z) 
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+ (c 1 {y{y + l) + {y + z) 2 ) 

-c 2 (2y + l)(y + z)) 
ln(-y) - ln(j/ 



+ 



+ 



1) 



g^y) 

(y + z) 2 
ln(z) — ln(z) 



2(y + z) 

(y + z? 



H(z,y) 



[y + zf 

+ {z^z). (19) 

4. Above formulas and the known NLO evolu- 
tion equations describing \ip dependence of the 
GPDs and meson DA give a complete basis for 
the description of a neutral vector meson electro- 
production with NLO accuracy. At leading order 
we reproduce the known result [9] , our results for 
the NLO hard-scattering amplitudes are new. 

At high energies, W 2 » Q 2 , the imaginary 
part of the amplitude dominates. Leading con- 
tribution to the NLO correction comes from the 
integration region £ <C |x| <C 1, simplifying the 
gluon hard-scattering amplitude in this limit we 
obtain the estimate 



M 



j*N^V L N 



(20) 
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Given the behavior of the gluon GPD at small 
x, F 9 (x,t;,t) <~ const, we see that NLO correc- 
tion is parametrically large, ~ ln(l/£), and nega- 
tive unless one chooses the value of the factoriza- 
tion scale sufficiently lower than the kinematic 
scale. For the asymptotic form of meson DA, 
4>y s (z) = Qzz, the last term in (20) changes the 
sign at fiF = g , for the DA with a more broad 
shape this happens at even lower values of fip. 

It is interesting to note that the value fj, F = 
Q 2 /e 2 is rather close to an estimate in the dipole 
approach [10] of a typical inverse dipole size 
[11,12], 1/r 2 ~ n 2 F = 0.15Q 2 , for vector me- 
son electroproduction in the HERA kinematic re- 
gion. We believe that a study of relationship be- 
tween the collinear factorization and the dipole 
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Figure 1. Factorization scale dependence of pre- 
dicted da L /dt\ t =o at Q 2 = 27 GeV 2 , W = 
110 GeV. The horizontal black line describes the 
data point. Solid curves assume [ir — [If, dashed 
curves - the BLM prescription. The curves a 
and c (b and d) are obtained with the use of 
MRST2001 (CTEQ6M). 



approach at the NLO level would be very impor- 
tant. But this question, as well as the resumma- 
tion of large at high energies, <~ (as ln(l/£)) n , 
contributions both to the hard-scattering ampli- 
tudes and to the evolution of GPDs, go beyond 
the scope of the present work. 
5. As an example of our results we com- 
pare a one point for the longitudinal cross sec- 
tion reported by ZEUS collaboration [13], e.g. 
da L /dt\ t =o = 17 ± 4 nb/GeV 2 at Q 2 = 27 GeV 2 , 
W = 110 GeV, with our predictions. On 
the Fig. 1 we plot the dependence of predicted 
do~L I dt\t=o(HF , Hr) on the factorization scale (jlf 
for two choices of the renormalization scale: /j,r = 
jip (the solid curves) and /in = Q/\/e, i.e. 
the BLM (Brodsky-Lepage-McKenzie) prescrip- 
tion (the dashed curves). The data point is de- 
scribed in this plot by the black horizontal line. 
In this numerical analysis we use two models of 
the NLO generalized parton distributions of [14] : 
first one based on MRST2001 forward distribu- 
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tion (the curves a and c) and second one - on 
CTEQ6M (the curves b and d). Moreover, we 
take the NLO strong coupling constant a s and 
the asymptotic meson DA 3 . Fig. 1 shows that the 
BLM prescription leads to smaller cross section 
and to much flatter behaviour of the cross-section 
on [if than for the choice pn = [if- We observe 
also a substantial uncertainty of our predictions 
due to the input parton GPDs. 

Since as mentioned above the NLO corrections 
are large it is instructive to study the relative 
magnitudes of different contributions to the am- 
plitude (10). This is done by assuming /j,r = fiF 
and for CTEQ6M GPD. In Fig. 2 we show plots of 
ImM and ReM as functions of \ip corresponding 
to a gluonic and a quark Born parts of the ampli- 
tude (10) (denoted as gB and qB, respectively), 
to the NLO part of the gluonic contribution T g 
(gN), to the NLO part of the quark contribution 
T q (qN), to the quark contribution T( + ) (q + ) and 
to the full amplitude (10) (denoted as full). All 
these separate contributions arc normalized by 
\ of (10) with all terms taken into account. Let 
us note that, as expected for the small x process, 
the imaginary part of the amplitude (10) domi- 
nates over its real part. We observe that the NLO 
corrections are mostly of opposite signs than the 
corresponding Born terms and they are big, con- 
sequently the final values of the amplitude (10) 
are the result of a strong cancellations between 
Born parts and NLO terms. Without account of 
NLO terms the predictions would be substantially 
above the data. These results are similar to ones 
obtained recently for T photoproduction [16]. 

In this paper we restrict our analysis to the 
leading twist and neglect completely the power 
suppressed, <~ l/Q, contributions. This is def- 
initely legitimate for sufficiently large Q. Most 
probably, at the HERA energies and Q 2 ~ 
20GeV 2 the higher twist corrections are still 
large, see e.g. [11,12] based on [17]. Nevertheless, 
we want to stress that our leading twist results ob- 
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Figure 2. Different contributions to ImM/\M\ 
and ReA4/\A4\ in dependence on fip (for fiR = 
Hf and CTEQ6M): gB - Born term of T g , qB - 
Born term of T q , gN - NLO terms of T g , qN - 
NLO terms of T q , q + - T (+) , full - all terms of 
(10) included. 



taincd with NLO hard-scattering amplitudes and 
NLO GPDs [14] ( which were adjusted to describe 
HERA deeply virtual Compton scattering data) 
are in qualitative agreement with the measured at 
HERA p meson electroproduction cross section. 
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